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Abstract 

We suggest that disordered framework aluminums and non-framework cations can 
create a disordered electrostatic potential in zeolites that can lead to a discrepancy be- 
tween diffusivities measured by microscopic and macroscopic experimental techniques. 
We calculate the value of the discrepancy and the characteristic time scale at which it 
occurs for several ionic and polarizable species diffusing in zeolites. For ionic species, 
a discrepancy is almost inevitable. For polarizable species, a significant discrepancy 
may occur in some zeolites only for long alkanes or large species such as benzene. 
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1 Introduction 

Disordered microporous materials, such as zeolites, microporous carbons, pillared clays, 
and superionic conductors, have a wide range of applications in fields such as catalysis, 
separation, remediation, and sensing. In all of these uses, the dynamical behavior of the 
sorbed molecules or ions can be a constraining factor. Due to the complicated guest-host 
interactions in these systems, our ability to predict the dynamical behavior on all time scales 
is still limited. Accurate prediction of diffusivities is, however, an important component in 
the design of practical processes using microporous materials. 

Among these materials, zeolites have been the focus of particularly intense research ef- 
forts due to their wide industrial use and unique crystal structures. Significant effort has 
been devoted to understanding diffusion in particular.! Despite this study, a fundamental 
question persists for some classes of zeolites: What causes the well-known discrepancy be- 
tween the diffusivities measured with macroscopic and microscopic techniques? Macroscopic 
methods, such as the zero-length column (ZLC) technique! and the Wicke-Kallenbach (WK) 
membrane method,! typically measure transport diffusivities under a concentration gradient. 
The time scales involved in the macroscopic experiments are long (~ Is). Microscopic tech- 
niques, such as NMR, measure self- diffusivities by following the motion of tagged molecules. 
In this case the system is kept under equilibrium conditions. The NMR technique samples 
motion on the Angstrom length scale and probes relatively short time scales (~ 10~ 12 -10~ 8 s) 
for the diffusivities typical in zeolites. Mesoscopic methods, such as NMR pulsed field gradi- 
ent (PFG)I and quasi-elastic neutron scattering (QENS),! also follow the microscopic motion 
of diffusing species. These methods are time-scale limited, able to access the 10~ n — 10~ 3 s 
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range. The diffusivities measured using macroscopic techniques can disagree from those 
measured by microscopic techniques by one or more orders of magnitude. For example, the 
diffusivity of benzene in NaX crystals at 403K is found to be 1 x 10 _10 m 2 s~ 1 by NMR, but 
4.5 x lO-^mV 1 by ZLC.i 

Several explanations for the discrepancy between macroscopic and microscopic diffu- 
sivities have been offered, including an extracrystalline resistance, surface barriers to mass 
transfer, and internal heat conduction effects.! To date, however, there is no general theory 
that can satisfactorily explain the discrepancy. Recent computer simulations have begun 
to shed light on the understanding of this long-standing puzzle. It appears that different 
types of molecular motions happen at different time scales. Molecular dynamics simula- 
tions by Maginn et al., for example, reveal that the time constant associated with molecular 
rearrangements of alkanes may be greater than the time constant for motion along a given 
channell Similarly, kinetic Monte Carlo simulation studies of benzene in zeolites by Auerbach 
et al. suggest that, due to the different experimental time scales, PFG NMR measures the 
intercage diffusivity, while the NMR relaxation method measures the intracage diffusivity! 

What these simulation studies suggest, and indeed what the discrepancy between mi- 
croscopic and macroscopic diffusivities suggests, is that the diffusivity is actually time de- 
pendent. That is, for some sorbates in some zeolites, there seems to be one value of the 
diffusivity, Dq, associated with motion on the Angstrom scale and another, smaller, value, 
-Dqo, associated with motions on longer spatial scales. Such a time-dependent diffusivity 
can arise from an unusual time dependence of the mean-square displacement of the diffusing 
species. In this paper, we provide a theory to calculate the diffusivity on different time and 
length scales. This theory gives not only the ratio D^/Do, but also the characteristic time 
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and length scale on which the measured diffusivity changes most dramatically. 

Although not generally appreciated, disorder within the zeolite framework can have a 
dramatic effect on the measured diffusivity. The aluminum framework species as well as the 
associated, non-framework cations create a disordered electrostatic environment for diffusing 
species if the aluminums or cations are randomly sited. The strength of the randomness is 
determined by the aluminum loading, which can vary from a 1:1 Si:Al ratio to a 3000:1 or 
higher Si:Al ratio. Higher levels of aluminum loading require higher densities of exchangeable 
cations.! Together, the framework aluminum and non-framework cations create a random 
environment and lead to a lower ratio of D^/Dq^ On the other hand, the disorder can either 
increase or decrease the microscopic diffusivity Dq. In the absence of any specific interactions 
between the diffusing species and the disordered ions, one would expect D to decrease with 
increasing disorder. For chemisorption-type interactions typical of nucleophilic species, Do 
may increase with an increasing density of cationic sites since the hopping distance between 
the sites decreases, yet the binding energy per site may remain constant.0 

In this paper, we model the long-distance features of the random electrostatic potential 
in the zeolite with a particularly simple, Gaussian form. We analyze this model to give an 
explicit form of the time-dependent, experimental diffusivity D(t)/D . With this result we 
can predict both the ratio of the macroscopic to the microscopic diffusivity and the time 
or length scale at which the difference between the two values becomes apparent. If this 
time scale is greater than that accessed by, say, the NMR technique, then disorder provides 
one possible mechanism for a discrepancy between microscopic and macroscopic diffusion 
measurements. 

To perform this analysis, we use a field theoretic formulation of the diffusion process. 
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While the theoretical treatment is involved, the physical mechanism considered is simply 
the slowing down of a diffusing species by the random electric fields within the zeolite. Our 
paper is organized as follows. In Section ^| we use renormalization group theory to derive the 
time-dependent diffusivity in terms of the disordered environment within the zeolite. The 
long-time renormalization group result is matched to a short-time perturbation theory to give 
an explicit prediction for the diffusivity. In Section |3] we describe Monte Carlo simulations 
that validate the time-dependent theory. Several forms of random disorder are studied, 
and some numerical issues are discussed. In Section £| we apply the theory to zeolites. The 
reader interested primarily in how disorder may explain the discrepancy between microscopic 
and macroscopic measurements of the diffusivity may wish to skip directly to this section. 
Diffusion of both ionic and polarizable species is considered in this section. Values of the 
discrepancy between microscopic and macroscopic measurements are predicted in several 
cases, and the characteristic time and length scales when the discrepancy becomes apparent 
are also given. For measurements of ionic diffusivities, discrepancies are almost inevitable. 
For polarizable species, discrepancies can be one order of magnitude or more for long alkanes 
or large species such as benzene. In Section |5|, we summarize our findings and make some 
final remarks. 



2 Time-Dependent Renormalization Group Theory 

Before specializing to the case of diffusion of ions or polarizable species at low densities 
in disordered zeolites, we first consider the more general problem of diffusion of any type of 
species in a disordered environment. In general, the disorder will slow down the diffusing 
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species, and we would like to predict the magnitude of this effect and the characteristic 
time scale at which it occurs. Since we will be interested only in the ratio of the macro- 
scopic diffusivity to the microscopic diffusivity, we can use a theory valid for long times and 
large distances. Field theory is just such a tool. Field theory also naturally allows us to 
consider the effect of the disordered environment. Alternatively, field theory can be viewed 
as one particularly effective way to solve for the Green function of this diffusion problem. 
Interactions between the diffusing species will be negligible at low densities, even for ionic 
speciesBIl! We, therefore, consider the single-species diffusion equation 



()a > {x -<) =DoV 2 G v (x,t)+f3D V- [G v {x,t)Vv(x)] (1) 



at 



Here v is the quenched random energy field caused by the disorder, G v is the Green function 
of a single diffusing species in a given realization of the disorder, D is the "bare," short-time 
diffusivity, and (3 = \/{hfeT) is the inverse temperature. The mean-square displacement is 
given by R%(t) = J d x |x| 2 G„(x, t). The mean-square displacement measured in an exper- 
iment is averaged over the statistics of the disorder: R 2 (t) = (R%(t)) = J d d y. |x| 2 G(x, t), 
where the Green function averaged over disorder is G(x, t) = (G v (x,i)). 

Because we are interested only in how the large-distance features of the random envi- 
ronment affect the diffusive motion, we can use a particularly simple, Gaussian form for the 
disordered energy field. Other forms turn out to be equivalent to this choice in a certain, 
technical sense. Since the force on a diffusing species is proportional to the gradient of the 
energy field, we can set the average value of the field to zero without loss of generality. 
And since we assume it is Gaussian, the only other quantity to specify is the energy-energy 



correlation function Xvvij) = ( v Wv(r)). 

The time-dependent field theory and model we use has been analyzed extensively in the 
context of reaction and diffusion processes in random media.0^1 Here we consider only the 
diffusion process. In this section, we derive a formula for the time-dependent self-diffusivity, 
which is defined as 

D(t) = (2) 
w 2dt K 1 

For the application to zeolites we will consider d = 3 dimensions, but here we develop the 
theory for a general spatial dimension. 

The derivation of the field theory is somewhat technical,ll3lll and so we provide a sum- 
mary of the procedure here. A master equation that describes the motion of the diffusing 
species is first written. This master equation is equivalent to eq |I[ The master equation is 
then rewritten in terms of creation and annihilation operators. The coherent state represen- 
tation is then used to rewrite the master equation in terms of a field theory. Finally, the 
replica trick is used to integrate out the effects of the disorder. After these exact transfor- 
mations, one arrives at the action S = S + Sj, with 

5 = EM'W dta a (x,t)[dt-D V 2 + 5(t)]a a (x,t) (3) 

a=l J J ° 

5 1 = tL lT- E / dt ^ / (2vr)^(ki + k 2 + k 3 + k 4 ) 

2 c ^L l JO Jk!k2k3k4 

a Q (ki, *i)a Q (k 2 , ti)a 7 (k 3 , t 2 )a 7 (k 4 , t 2 ) 

kx • (k! + k 2 )k 3 • (k! + k 2 )x OT (|k 1 + k 2 |) (4) 
The Fourier transform of a function is defined by /(k) = / c/ d x /(x) exp(zk • x), and the 
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notation J* k is a short hand for / d d ~k / (2ft) d . The number of replicas is N, and the iV — > 
limit is implied. The action is used to calculate averages. The Green function, for example, 
is given by 

1 N 

G(x,t) = lim-^(a 4 (x,t)a t (0,0)) (5) 

i=i 

where the average is with respect to the weight exp(— S). The free theory term So represents 
simple diffusion. The interaction term Si comes from averaging the diffusion process over the 
random energy field, which is again assumed to be Gaussian, with zero mean and correlation 
function \ vv . The form of \ vv is left general throughout the calculation. 

Using field theoretic renormalization group theoryS we let a — > Z a 2 a, a — ► Z£a, D — ► 
Z~ 1 ZdD, and j3 2 — > /3 2 Zp Zj. Note that the dimensions of P 2 Xw are t ne same as that of 
k~ d , where k has the dimensions of wave number. We have chosen f3 as dimensionless in 
this renormalization scheme. The Z factors will be determined by a subtraction scheme. 
Following a one- loop calculation of the renormalized vertex functions in the long-time limit, 
we find 

Z a = 1 

Z D = [ A dk k d ~ l x vv (k) 

Zi = l- 2 -^ J\kk d - l x m {k) (6) 

where K d = S d /{27i) d , S d = 27r d / 2 /[r(d/2)] is the surface area in d dimensions, A is a cut-off, 
and \i = A/e l with / the flow parameter. No contribution renormalizes the time derivative 
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or delta function term in S . The associated renormalization-group functions are defined by 



1% = V 



dfx 



, i = a,D,I 



(7) 



where the subscript on the derivative indicates fixed bare variables. We find 



7a = 



7d = (3 2 K d n d x vv (fi)/d 



7/ = 2p 2 K d ii d x m (v)/d 



and 



Pi = fJ> 



8f3 2 



d/j, 



= (7/ - 2 lD )f3 2 = 



(9) 



The dynamical exponent is given by 



z = 2- la + lD = 2 + [3 2 K d fi d Xw(v)/d 



(10) 



The renormalized time flows according to 



t(l*) = texp 



- [ l z(l)dl 
Jo 



(11) 



The matching time, to, is usually on the order of 1/(A 2 _D ), which is within the range of 
validity of both renormalization group scaling and mean field theory. The Callan-Symanzik 
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equation then reads 



(12) 



Solving this differential equation, we find 



D(l*) = D(0)exp 



(3 2 K d fi* (A\\ /A 

CM I — I X 



d Jo 



(13) 



An explicit form for D(t) follows from eq [13] after relating the flow parameter I* to time by 
eqs 



10 and 11 



D(t) = D exp 
= D exp 



(3 2 K d f'*W./A\^ /A 
di I — I x 



d Jo 

(3 2 K d rA 



d J\/e l 



1 I /\VV \ 1 



Jkk d ~ l x vv {k) 



(14) 



In the limit of t — ► oo, which implies I* — > oo, the above equation yields 



D oo = D ex P [-p 2 X vv(0)/d} 



(15) 



which is in agreement with the long-time renormalization group result found in Refs |T0 and 



17|. Finally, eq 11 can be rewritten as 



t = t e 2l ~D /D(t) 



(16) 



The renormalization result depends on the matching time to appearing in eq 11. We 



determine to by equating perturbation theory to the renormalization group result for weak 



10 



disorder, or equivalently, (3 — > 0. The mean square displacement is given by R 2 (t) = 2dDt, 
so the diffusivity is related to the Green function by 



D(t) 



Jd d xx 2 G(x,t) V£G(k,t)| 



fc=0 



2dt 



2dt 



(17) 



For weak disorder, perturbation theory gives 



D(t) = D - / A dk k d - 3 x vv (k)(-3 + x + 3e~ x + 2xe~ x ) + 0(/? 4 , t 3 ) (18) 

dt Jo 



where x = D k 2 t. In the limit (3 — > we match perturbation theory and the renormalization 
group theory at a characteristic time, t c : 



D(t c ) = (D + D 00 )/2 = D c 



2u 



(19) 



We find t c by equating eqs ITBI and FBI to 0(/5 



Xw(0) = —± I dk k d - 3 Xw(k)(-3 + x c + 3e~ x ° + 2x c e~ Xc ) 

Dntr JO 



D t c Jo 



(20) 



where x c = D^k t c . Note that eq 20 shows that t c is a function of A, Dq, and the form of 



Xw, but not the strength of the disorder. We similarly find l c by equating eqs [14] and [19] to 

o(p 2 y. 

Xvv (0) = 2K d dl xw (^) (21) 



We finally obtain to at small (3 from the relation, to = t c e . We assume that to is only 
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weakly dependent on the strength of the disorder and use the same t for all values of (3 2 . 



3 Monte Carlo Simulation Results 

Before we apply the time-dependent theory to diffusion in disordered microporous ma- 
terials, we first test it against computer simulations. We do this by generating numerical 
solutions to eq [L| averaging them over the statistics of the disorder, and then comparing 



the results to eq 14, The computer simulations are done on a lattice in two dimensions and 



follow closely the strategy of Ref |18[ Using eqs (3-10) from Ref [18], we create a Gaussian 
random energy field, prescribed by v(k), on a lattice. We determine the energy field in real 
space by fast Fourier transform. To obtain reasonable statistics for the diffusivity, we follow 
the motion of 10000 independent particles on a lattice of size L. Periodic boundary condi- 
tions are enforced. We use lattices of size L = 512, 1024, and 2048 in order to identify and 
avoid finite size effects. The initial particle positions are chosen at random. Finally, D(t) is 
measured by averaging R 2 (t) as a function of t over a moving window of 200 Monte Carlo 
steps and evaluating the slope. To test the robustness of the theory, we examined three 
forms of Xvv- Xw = 7 ^p[~k 2 / (2k 2 )}, Xw = 7exp(-|A;|/A; c ), and Xw = lki/{k 2 + k 2 ) 2 . 

The simulation results with L = 2048 and k c = 0.1 are shown in Figure [l]-[|for each form 



of Xvv Significantly, the infinite time result for the diffusivity, eq [L^, appears exact to within 
the statistics of our simulation. More importantly for the present purpose, there is a good 
agreement between the renormalization group result, eq O, and the simulation data over 



the entire time range. Indeed, to appears to be only weakly dependent on the strength of 
disorder. Recall that to was determined at weak disorder, and we assumed it to be constant 
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at strong disorder. The observed constancy of t validates our matching procedure. 

Interestingly, two types of finite size effects can occur in these simulations. First, the 
lattice must be large enough to capture the length scale at which the disorder begins to slow 
down the diffusing species. If the lattice is not large enough, the diffusivity will reach the 
proper asymptotic, long-time value. Second, the lattice must be fine enough so as to capture 
all of the features inherent in the chosen ^ TO (x). If the characteristic length of Xm>( x ) is 
finer than the lattice spacing, then some of the effects of the random energy field will be 
lost. Since the diffusion constant decreases with increasing disorder, a lattice that is not fine 
enough and, therefore, neglects some of the effects of the disorder, will predict a diffusivity 
that is too large. We, indeed, observe this effect for simulations with k c = 1.0. In this case, 
the characteristic size r c = 2n/k c is approaching the lattice spacing of unity that we use. 

Note that we observe a peak of D(t)/Do in the simulations for small D^k^t. This peak 
is more evident for small values of k c . In Figure | we plot the simulation results together 
with the renormalization group and short-time perturbation theory results for k c = 0.02. It 
is clear that perturbation theory describes the short-time peak well. This is to be expected 
since the perturbation result is accurate to 0(t 2 ) for all values of (3 2 . In contrast, the 
renormalization group approach fails to predict this interesting short-time behavior, since 
the renormalization procedure is carried out in the asymptotic, long-time limit. This peak 
in the diffusivity is due to the rapid relaxation of the particles away from the initial, non- 
equilibrium distribution. Consequently, this peak would not be expected to be observed in 
equilibrium techniques such as NMR, PFG NMR, or QENS. 
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4 Time-Dependent Diffusion in Zeolites 

Within a zeolite crystal there is a random potential field, mainly due to the random 
electrostatic potential created by the charged framework Al _ and non-framework cation 
sites. The lattice defects, cations, and associated aluminum atoms create strong disorder, 
increase the hopping activation energy, inhibit the mobility of sorbates, and slow down the 
diffusion process. We list the ionic self-diffusivities measured in same-type of zeolites with 
different aluminum loadings in Table [I]. To reduce the effects of ionic hydration, data taken 
for methanol or ethanol solvents are exhibited. As can be seen, Table [l] confirms a reduction 
diffusivity caused by aluminum-induced disorder. That is, at higher aluminum loadings the 
activation energy increases and the diffusivity decreases. 

We assume that the spatial arrangement of the charges in zeolites is fixed during for- 
mation conditions, i.e. that the disorder is quenched.0 The Debye-Huckel theory^ for the 
density-density correlation function of the charged impurities then gives 

X PP (k) = (22) 

with 7 = |e~| 2 (p + + and k is an inverse characteristic length of the charged disorder, 
defined by k 2 = j3 F \e~\ 2 (p + + p_)/e F . Here e~ is the electron charge, p + and p_ are the 
positive and negative charge number densities, respectively, ep is the dielectric constant 
under formation conditions, and (3-p = I/^b^f) with Tp the formation temperature. Under 
typical aqueous synthesis conditions of pH = 12 and Tp = 398 K, the dielectric constant cf 
is approximately €f = 5Oeo-0 The correlation function of the electrostatic potential due to 
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the charged impurities is: 



*#(*) = (23) 



where now the dielectric constant should be the value in the zeolitic solid, i.e. e = 4.2eo-0 

We first consider the diffusion of an ion of charge q in such a zeolite. The energy-energy 
correlation function due to the charged disorder is 

xZ(k) = q 2 X*m (24) 

We consider diffusion of a singly charged species, q = \e~\, at room temperature. We take 
zeolite NaY with px = 1-27 x 10 28 m -3 and pai = 0.3 px as our example, where px is the 
density of T (=Si or Al) atoms in the zeolite. We are interested in the magnitude of the 
time-dependent diffusivity as well as the characteristic time, t c . It appears at first sight that 
our theory requires no cutoff. Directly solving eq ^ and ^D], we find D^/Dq = e -634 , and 



D K 2 t c = 8.8. The term /3 2 x{fy/d in the exponential of eq 15 looks like an activation energy 
introduced by disorder, but the numerical result, (3E = 634, is significantly larger than that 
measured in experiments.il 

It is, however, clear that Debye-Hiickel theory is not quite right for short distances or 
large wave vectors. In particular, the characteristic length of the charged disorder, 2tt/k, 
is typically on the order of several Angstroms. Motion on the Angstroms length scale is 
measured in microscopic experiments such as NMR and is well described by Newton's equa- 
tion, e.g. by molecular dynamics simulations. Statistical field theory is applicable only for 
behavior at much longer time or length scales, and so Dq is an input parameter, determined 
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from a simulation or a short-time experiment. Since this D already captures short time 
and distance properties, a cutoff is required to exclude these effects. This cutoff will be a 
function of k. Given the measured diffusivity and the experimental time scale, t exp , however, 
the cutoff can be calculated by 

A-exp 27r/r CX p 

r 2 cxp = 2D t exp (25) 



Physically this means that we should use field theory only to consider fluctuations with 
wavelengths larger than r exp . 

Many factors, such as charge neutrality, Loewenstein's rule, strong cation hydration, 
ion-sieving effects, and the existence of several different energetically favorable cation sites, 
all complicate the cationic diffusion process and make ion diffusion data hard to organize.0 In 
Table |2], we illustrate the effect of disorder for a few typical diffusivities. Since the available 
techniques for measuring cation diffusion in zeolites, such as tracer exchange or counter 
diffusion exchange, are at the time scale of minutes, we use t exp = 100 s. From Table [| 
we see that a smaller D leads to a larger cutoff and a larger correction to the diffusivity. 
For example, if D = 1 x 10 -10 m 2 s -1 we see D^/Dq = 1 and no correction, while with 
Dq = 1 x 10 _17 m 2 s _1 we see a significant correction of Doo/Dq = 0.002. For all these 
cases, the cutoff A exp obtained is much smaller than k, and so the field-theoretic approach 
is appropriate. By assuming A < k, we can further simplify the energy-energy correlation 
function as 

2 

Xw(k) = ^; with rj = |-J (26) 
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From this form of x, we find 



Doo = D Q exp[-/3 2 K d r]A exp /d} 



(27) 



Solving for t c from eq 20, we find 



Al p D t c = 37.3 



(28) 



and from eq 21, we find 



to = t c /4 



(29) 



From eq [2^ we see immediately that t c oc t exp , which explains why we would have the 
same value of t c for the different diffusivities in Table 0. When the diffusivity is dramatically 
reduced, the characteristic time is better defined as the time at which the measured diffusivity 
is very close to Doo, i.e. D(t 90 %) = lAD^. We have calculated t$ % numerically from eqs [H 
and [16] and list the values in Table [|. As expected, t 90 % increases with decreasing D^/Dq. 
Interestingly, we find the mean-square distance that ions diffuse within this time period is 



constant: 



Ax 90 % = [R 2 (t 90 %)/3] 



1/2 



P 2 V 



3vr ln(l.l) 



.1/2 



exp _ 



(30) 



We now turn to consider the case of diffusion of polarizable species in zeolites. The 
energy of a polarizable species in an electric field is given by v (x) = — |a|V0(x)| 2 , where 
a is the molecular polarizability, and cf) is the electric potential. Using the Gaussian ap- 



proximation for 0, we find the energy-energy correlation function of polarizable species to 
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be 

xL° larizaWc (k) = \o? / p [(p - k) ■ p] 2 x^(p - k)**(p) (31) 

A cut-off is necessary to ensure the convergence of this integral. As before, this cutoff must 
be smaller than k if the field theory is to be applicable. Furthermore, the length scale should 
also be on the order of or larger than a molecular size and roughly equal to the displacement 
of molecules on the microscopic experimental time scale. We simply set A = 27r/(5A). 

We list results of the calculation for the four species H 2 , CH 4 , C 4 H 10 , and C 6 H 6 in Table 
H]. For H 2 , the transport and self-diffusion coefficients were measured simultaneously in Ref |5] 
by QENS. No discrepancy was found, and the diffusivities given by the two methods are in 
good agreement with each other, D = 3 x 10~ 9 m 2 s -1 at T = 100 K in NaX. Methane 
serves as another example where no discrepancies have been reported. The diffusivity of 
methane in NaX measured by QENS is 3.2 x 10 _9 m 2 s -1 .il Our calculations support the 
lack of discrepancy in these cases. Although the characteristic time scales lie just within 
the range of QENS experiments, the calculated D^/Dq for H 2 and CH4 is 0.72 and 0.39 
respectively. Such a ratio would be difficult to observe, given the small characteristic times. 

Numerous experiments have reported that there are two orders of magnitude of dis- 
crepancy between the microscopic and macroscopic diffusion of benzene in NaX or NaY.i 
Computer simulations by Auerbach et al. of benzene diffusion in Na-X or NaY suggest that 
the longer time scale PFG-NMR technique measures intercage diffusion, while the short time 
scale NMR spin-lattice relaxation technique measures intracage diffusion. These simu- 
lations also suggest that the cage-to-cage length, 11 A, is a reasonable approximation for the 
effective hopping length. Our cutoff should be smaller than this value to capture the effects 
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of disorder on this length scale. We use the microscopic D = 1.9 x 10 _10 m 2 s _1 measured 
by the QENS method.il With our chosen cutoff, our theory predicts t 90 % = 4.4 x 10~ 8 s 
and Doc/ Do = 0.08. The large characteristic time is primarily due to the smaller diffusivity 
for benzene compared to hydrogen or methane. For this species, the greater discrepancy in 
diffusivity is primarily due to the greater polarizability of benzene compared to hydrogen 
or methane. For this species, we may expect to see the time-dependence of the diffusiv- 
ity in QENS experiments, or we may expect to see discrepancies between microscopic and 
macroscopic measurement s.S Consistent with the idea that increased disorder slows down the 
diffusion species, benzene diffusivities observed in simulations also decrease with increased 
aluminum loading.0 For butane, the calculated D^/Dq is small. Since the characteristic 
time £90% lies within the range of QENS, we may be able to observe a discrepancy in this 
case as well. 

Instead of imposing a length-scale cutoff for diffusion of the polarizable species, we may 
impose a time-scale cutoff, as we did in Table ^. We choose t exp = 10~ 9 s, since this is at 
the lower end of the time scale accessible to QENS and near the upper end of the time scale 



accessible to molecular dynamics. We use eq [25| to derive the value of A cxp to use in our 
calculations. The results are listed in Table |3|. With this type of cutoff, we find a discrep- 
ancy between microscopic and macroscopic measurements only for benzene. Comparison of 
the results in Tables [| and shows that differences between microscopic and macroscopic 
measurements depend sensitively on the details of the time or length scale resolved by the 
microscopic experimental technique. Experimental diffusivities may be less sensitive to the 
cutoff procedure than these tables would suggest, since the Gaussian approximation that we 
have used has lead to an over- sensitivity to short-distance features in the case of polarizable 
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species. 

On the other hand, diffusion data for many species in zeolite A are in satisfactory 
agreementJl This seems at odds with our theory, since in zeolite A, the Si:Al ratio is as 
high as unity. However, a closer examination of the structure of zeolite A reveals a perfect 
crystal structure. The aluminum ions obey Loewenstein's rule and are arranged in a highly 
symmetrical way, such that each unit cell consists of 96 AIO4 and 96 SiC>4 tetrahedra with 
each apical oxygen atom shared between one aluminum atom and one silicon atom. Zeolite 
A is not disordered, and so the movement of molecules within each unit cell rapidly reaches 
diffusive behavior. We, therefore, do not expect to observe any long time corrections to 
diffusivities in zeolite A. Although the Si:Al ratio is close to unity in zeolite X as well, 
the unit cell is significantly larger and the electrostatic potential is substantially disordered 
compared to zeolite A, as judged by the presence of the strong electric field gradients that 
make X-type zeolites suitable for selective separation of air.! It is this disorder that can lead 
to a time-dependent diffusivity. 

We can summarize our results by noting that the characteristic time for the slowing 
down of the diffusing species due to disorder is given by t c = £ 2 / (2D ), where £ is a charac- 
teristic length associated with the disorder. Smaller values of D , therefore, lead to larger 
characteristic times and a greater chance that a discrepancy between microscopic and macro- 
scopic measurements will occur. We have considered only two contributions to the force on 
a diffusing species, the Coulombic force on ions and the induced-dipole force on polarizable 
species. Other interactions between the diffusing species and the zeolite framework and 
non-framework atoms, such as ion-dipole, ion-quadrapole, Lennard- Jones, and so on, also 
contribute to the random energy fields. These additional terms would also slow down the dif- 
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fusing species, although the characteristic times of most of these interactions are too short to 
cause an observable discrepancy between microscopic and macroscopic measurements. Our 
calculations for diffusion of ions should be fairly robust, since in this case there is especially 
good reason to believe that the disorder is well-approximated by a Gaussian.0 For polar- 
izable species, our calculated corrections to the short-time diffusivity are sensitive to the 
numerical value of the cutoff. Equivalently, the theory is sensitive to the assumptions about 
the behavior of the dielectric constant at short distances. We have also assumed that the time 
scales associated with configurational arrangement of the diffusing, polarizable species are 
much shorter than the characteristic diffusion time induced by the disorder. Our predictions 
for the ionic case, therefore, are more quantitatively reliable than those for the polarizable 
case. Regardless of the details, that disorder induces a time-dependent diffusivity, as in eq 
refeq2, for an ionic or polarizable species is a general result. 

Disorder on longer mesoscopic length scales can also occur in zeolites in the form of 
faulting or growth-induced domain-mismatch. The length scales associated with such do- 
mains can be on the order of £ = 10 nm in the ZSM-fS and Linde Type L0 zeolites. It 
would be interesting to see if a reduction in the diffusivity occurring at t c = (, 2 /(2D) is 
observable experimentally for these materials. 

5 Conclusion 

We have derived a time-dependent theory of diffusivity in disordered microporous me- 
dia. As judged by computer simulations, the theory is reasonably accurate for Gaussian-type 
disorder over the entire time range. Applying the theory to zeolites, we offer disorder as one 
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possible mechanism for the existence of a discrepancy between microscopic and macroscopic 
measurements of diffusion in zeolites. For diffusion of ionic species, we predict that a dis- 
crepancy almost certainly exists. For polarizable species such as H 2 and CH 4 we predict 
no observable discrepancy, and none has been observed experimentally. Our study of ben- 
zene and n-butane suggests that disorder may cause a discrepancy for these molecules. The 
characteristic time and magnitude of the discrepancy are consistent with microscopic and 
macroscopic experimental observations. We explain the lack of a discrepancy for zeolite A 
as resulting from a lack of disorder in this highly-symmetric, perfectly-crystalline material. 

Our theory can only predict the ratio D(t)/D . Determination of the short-time dif- 
fusivity, Do, is best done by a molecular dynamics simulation or a microscopic experiment. 
In this respect, molecular dynamics simulation is a particularly useful, complementary tool 
to the field theory presented here. Diffusivities derived by molecular dynamics are often 
comparable with the NMR, PFG NMR, and QENS measurements.il Moreover, simulation 
provides molecular details beyond the scope of an experiment or a field theory. Such simu- 
lations should, for example, allow one to construct a more exact form of x vv (k) to be used 
in our theory. 
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Table 1: Measurements of self diffusivities and activation energies in zeolite Y with 
different Al content.il 



Zeolite 


Si:Al 


Pai/ptO 


Solvent 


E(kJ/mol) 


D (m 2 s 




Sr-Y 


139:53 


0.276 


Methanol 


91.6 


5.2 x 10 


-22 


Sr-Y 


125:67 


0.349 


Methanol 


118 


1.4 x 10 


-25 


Ba-Y 


139:53 


0.276 


Methanol 


67.8 


3.3 x 10 


-20 


Ba-Y 


125:67 


0.349 


Methanol 


84.1 


1.5 x 10 


-21 


Ba-Y 


139:53 


0.276 


Ethanol 


91.6 


3.3 x 10" 


-24 


Ba-Y 


125:67 


0.349 


Ethanol 


120 


3.1 x 10" 


-25 



a pT is the density of framework T (=Si or Al) atoms. 



Table 2: Ionic diffusivities in NaY with pai = 0.3 px and T = 298 K. 



D (m 2 s 


- 1 ) 


^exp (s) 


A cxp (m x ) 


Doc/Do 


^90% (s) 


Ax 90 % (m) 


1 x 10- 


10 


100 


4.4 x 10 4 


1 






1 x 10- 


14 


100 


4.4 x 10 6 


0.82 


230 


2.03 x 10~ 6 


1 x 10" 


15 


100 


1.4 x 10 7 


0.53 


3.5 x 10 s 


2.03 x 10 _B 


1 x 10- 


16 


100 


4.4 x 10 7 


0.14 


1.4 x 10 5 


2.03 x 10~ 6 


1 x 10- 


17 


100 


1.4 x 10 8 


0.002 


1.0 x 10 8 


2.03 x 10~ 6 


1 x 10- 


18 


100 


4.4 x 10 8 


e -20 


8.6 x 10 14 


2.03 x 10~ 6 



Table 3: 



Characteristic time, length, and D^/Dq for polarizable species with A^ 27r/(5A).g 



Species 


Zeolite 


PAl/PT 


T(K) 


D (m 2 s 


- 1 ) 


Ref a/(4vre ) (m 3 )EI 


Doo/Do 


^90% ( s ) 




Ax 90 %(A) 


H 2 


NaX 


0.5 


100 


3 x 10- 


9 




I 


0.8 x 10~ su 


0.72 


7.1 x 10- 


11 


5.8 


CH 4 


NaY 


0.3 


150 


3.2 x 10 


-9 


23 


2.6 x 10^ u 


0.39 


3.0 x 10- 


10 


9.1 


C4H10 


NaX 


0.5 


298 


2.5 x 10 


-9 







8.2 x 10~ 3U 


0.02 


2.0 x 10" 


-8 


15 


C6H6 


NaX 


0.5 


458 


1.9 x 10" 


-lu 




25 


10 x 10- 3U 


0.08 


4.5 x 10" 


-8 


13 



"If a short-distance dielectric constant of e = 2.0eo were used, a cutoff of A = 27r/(9.5A) would 
lead to similar values of Doo/Do and characteristic times increased by roughly a factor of three. 



Table 4: 



Characteristic time, length, and D^/Do for polarizable species with t exp = 10 9 s|] 



Species 


Zeolite 


Pax/ Pi 


T(K) 


Do (mV 1 ) 


Ref a/(4^e ) (m 3 )EI 


Doo/Dq t 90 % (s) Ax 90 %(A) 


H 2 


NaX 


0.5 


100 


3 x 10" 9 






0.8 x 10~ 3U 


1 


CH 4 


NaY 


0.3 


150 


3.2 x 10" 9 


1 


23 


2.6 x 10" 3U 


1 


C4H10 


NaX 


0.5 


298 


2.5 x 10~ 9 






8.2 x 10~ 3U 


0.999 


C6H6 


NaX 


0.5 


458 


1.9 x 10~ 1U 




25 


10 x 10- 3U 


0.37 7.5 x 10~ 9 11 



a If a short-distance dielectric constant of e = 2.0eo were used, a i Gxp = 3.5 x 10 9 s would lead to 
similar value of Doo/Dq and a characteristic time increased to 1.9 x 10 _8 s for benzene 
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ln(D k t) 
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Figure 1: Comparison of the theoretical results (dashed) with simulation data (solid) for 
Xw(k) = 7exp[— k 2 /(2k 2 )}, with k c = 0.1. Three different values of 7 were used correspond- 



ing to Doc/Do = 0.8, 0.5, and 0.1, as predicted by eq [15 



0.9 



0.7 



9 
a 



0.5 



0.3 




0.1 



-0.5 



1.5 



3.5 

ln(D k" 



5.5 



7.5 



Figure 2: Comparison of the theoretical results (dashed) with simulation data (solid) for 
Xw(k) = 7exp(— \k\/k c ), with k c = 0.1. Three different values of 7 were used corresponding 
to Doo/Dq = 0.8, 0.5, and 0.1, as predicted by eq |T5|. 
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-0.5 1.5 3.5 5.5 7.5 

ln(D k e t) 



Figure 3: Comparison of the theoretical results (dashed) with simulation data (solid) for 
Xw(k) = •yk^/(k 2 + k 2 ) 2 , with k c = 0.1. Three different values of 7 were used corresponding 
to Doo/Dq = 0.8, 0.5, and 0.1, as predicted by eq |T5|. 




ln(D k/t) 



Figure 4: Comparison of the renormalization group (dashed) results with simulation data 
(solid) and short-time perturbation theory (long-dashed) for Xw{k) = 7exp[— k 2 /(2k 2 )], with 
k c = 0.02 and D^/Dq = 0.1. 
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